Deterministic entanglement of two trapped ionstl 

Q. A. Turchettejll C. S. Wood, B. E. King, C. J. Myatt, D. Leibfried,! W. M. Itano, C. Monroe, and D. J. Wineland 

Time and Frequency Division, National Institute of Standards and Technology, Boulder, CO 80303 

(February 1, 2008) 



00 

On 
G 

1—5 

m 



> 

(N 

o 

\o 
o 
oo 

a,: 

c ■ 

3 : 
cr 



X 



We have prepared the internal states of two trapped ions in 
both the Bell-like singlet and triplet entangled states. In con- 
trast to all other experiments with entangled states of either 
massive particles or photons, we do this in a deterministic 
fashion, producing entangled states on demand without selec- 
tion. The deterministic production of entangled states is a 
crucial prerequisite for large-scale quantum computation. 



Since the seminal discussions of Einstein, Podolsky, 
and Rosen, two-particle quantum entanglement has been 
used to magnify and confirm the peculiarities of quantum 
mechanics Jlj] . More recently, quantum entanglement has 
been shown to be not purely of pedagogical interest, but 
also relevant to computation [^) , information transfer || , 
cryptography Q and spectroscopy Quantum com- 

putation (QC) exploits the inherent parallelism of quan- 
tum superposition and entanglement to perform certain 
tasks more efficiently than can be achieved classically . 

Relatively few physical systems are able to approach 
the severe requirements of QC: controllable coherent 
interaction between the quantum information carriers 
(quantum bits or qubits), isolation from the environment, 
and high-efficiency interrogation of individual qubits. 
Cirac and Zoller have proposed a scalable scheme uti- 
lizing trapped ions for QC ||. In it, the qubits are two 
internal states of an ion; entanglement and computation 
are achieved by quantum logic operations on pairs of ions 
involving shared quantized motion. Previously, quantum 
logic operations were demonstrated between a single ion's 
motion and its spin the requirements of QC have been 
explored experimentally in related cavity QED systems 
p0| | . In this Letter, we use conditional quantum logic 
transformations to entangle and manipulate the qubits 
of two trapped ions. 

Previous experiments have studied entangled states of 
photons ||n] , |l2i and of massive particles These 
experiments rely in some way on random processes, ei- 
ther in creation of the entanglement, as in photon cas- 
cades |OL photon down-conversion |l^] and proton scat- 
tering W3|, or in the random arrival times of atoms in 
a cavity [Hj. Recent results in NMR of bulk samples 
have shown entanglement of particle spins |fl5| , |l6[ but be- 
cause pseudo-pure states are selected through averaging 
over a thermal distribution, the signal is exponentially 
degraded as the number of qubits is increased. All the 
above processes are selectable but are not deterministic 
generators of entanglement. By deterministic, we mean 



that a known and controllable quantum state of (all of) a 
given set of particles is generated at a specified time [p.7| . 
Deterministic entanglement coupled with the ability to 
store entangled states for future use is crucial for the re- 
alization of large-scale quantum computation. Ion-trap 
QC has no fundamental scaling limits; moreover, even 
the simple two-ion manipulations described here can, in 
principle, be incorporated into large-scale computing, ci- 
ther by coupling two- ion subsystems via cavities Jl8|| , or 
by using accumulators || . 

In this Letter, we describe the deterministic generation 
of a state which under ideal conditions is given by 
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where | J.) and | \) refer to internal electronic states of 
each ion (in the usual spin- 1/2 analogy) and cf) is a con- 
trollable phase-factor. For <j> = or ir, \ip e ((f)) is a good 
approximation to the usual Bell singlet (— ) or triplet 
(+) state = [| |t) T I U)]/V2 since |(V s |<M0))| 2 = 
\{iPb\Mk))\ 2 = 0.98 jig. The fidelity of our experi- 
mentally generated state described by density matrix p 
is (^ e ( 7 r,O)|p ± |'0e(7r,O)) » (V'||P ± I^|) « 0.70, so that 
for all practical purposes, we can consider p ± to be an 
approximation to the Bell states. We describe a novel 
means of differentially addressing each ion to generate 
the entanglement and a state-sensitive detection process 
to characterize it. 

The apparatus is described in Ref. |pC|j , We confine 
9 Bc + ions in an elliptical rf Paul trap (major axis « 
525/im, aspect ratio 3:2) with a potential applied between 
ring and end-caps of Vb cos51t< + C^o with U,t/2tt ss 238 
MHz, Vq « 520 V. The trap is typically operated over the 
range 12 V < Uo < 17 V leading to secular frequencies of 
K,w y ,w z )/2jr = (7.3,16,12.6) to (8.2,17.2,10.1) MHz. 
The ion-ion spacing (along x) is I ~ 2/xm. 

The relevant level structure of 9 Be + is shown in Fig. 
|l|a. The qubit states are the 2s 2 Si/2 \F — 2,mp — 

2) = | i) and 2s 2 S 1/2 \F = l,m F = 1) = | T) states. 
Laser beams Dl and D2 provide Doppler precooling and 
beam D3 prevents optical pumping to the \F — 2,mF — 
1) state. The cycling | j) -> 2p 2 P 3/2 \F = 3,m F = 

3) transition driven by the c + -polarized D2 laser beam 
allows us to differentiate | f) from | |) in a single ion with 
«90% detection efficiency by observing the fluorescence 
rate. 

Transitions | l)\n) «-» | t)l n ') (where n,n' are vibra- 
tional quantum numbers) are driven by stimulated Ra- 
man processes from pairs of laser beams in one of two 
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geometries. Two types of transitions are driven: the "car- 
rier" with n' = n, and the red motional sideband (rsb) 
with n! = n— 1 m. With reference to Fig. 0a, the pair of 
Raman beams Rl _L R2 has difference- wavevector Sk \\ x 
and is used for sideband cooling (to prepare | ||)|0)), 
driving the i-rsb, and to drive the "i-carrier" . Beam 
pair R2 || R3 has Sk « and is not sensitive to motion; 
this pair drives the "co-propagating carrier" transition. 

Two trapped ions aligned along x have two modes 
of motion along x: the center-of-mass (COM) mode 
at frequency lo x and the stretch mode, at frequency 
u! stI — \fiuj x in which the two ions move in opposite di- 
rections. We sideband-cool both of these modes to near 
the ground state, but use the stretch mode on transitions 
which involve the motion since it is colder (99% probabil- 
ity of \n = 0)) than the COM and heats at a significantly 
reduced rate [^0) . The relevant two-ion qubit level struc- 
ture dressed by the quantized harmonic stretch motion 
is shown in Fig. 0b (we leave out the COM for clarity). 
In general, all four Rabi rates f2j±, i € {1, 2} connecting 
the levels are different and depend on n. Fig. 0b shows 
the states coupled on the rsb with Rabi frequencies (in 
the Lamb-Dicke limit) 
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where vj = 77/v2a/3 is the stretch- mode two-ion Lamb- 
Dicke parameter (with single-ion r\ w 0.23 for lo x /2tt « 
8 MHz) and fij is the carrier Rabi frequency of ion i 
||. On the carrier, the ions are not coupled and the 
time evolution is simply that of independent coherent 
Rabi oscillations with Rabi frequencies Qj. On the co- 
propagating carrier, fix — fl 2 = O c . 

In the Cirac-Zoller scheme, each of an array of tightly 
focused laser beams illuminates one and only one ion for 
individual state preparation. Here we pursue an alter- 
native technique, based not on f2, — ■+ for all but one 
ion, but simply on fii ^ Cl 2 . Differential Rabi frequen- 
cies can be used conveniently for individual addressing 
on the i-carrier: for example, if fii = 2Q 2l then ion 1 
can be driven for a time f2]i = tt (27r-pulse, no spin-flip) 
while ion 2 is driven for a 7r-pulse resulting in a spin-flip. 

Our technique for differential addressing is to control 
the ion micro-motion. To a good approximation, we can 
write 
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where Jo is the zero-order Bessel function and £j is the 
amplitude of micro-motion (along x) associated with ion 
i, proportional to the ion's mean displacement from trap 
center. The micro-motion is controlled by applying a 
static electric field to push the ions |2^] along x, moving 
ion 2 (ion 1) away from (toward) the rf null position, 
inducing a smaller (larger) Rabi frequency. The range 
of Rabi frequencies explored experimentally is shown in 
Fig. §i. 

We determine f2j 2 by observing the Rabi oscillations 
of the ions driven on the ii-carrier. An example with 



S7i = 2£l 2 is shown in Fig. gb. We detect a fluorescence 
signal S(t) = 2P U + (1 + a)P iT + (1 - a)P n where P M = 
\(ip(t)\kl)\ 2 , k,l € {T,l}, i>(t) is the state at time t and 
\a\ <C 1 describes a small differential detection efheiency 
due to the induced differential micro-motion. Driving on 
the ai-carrier for time t starting from | j|)|0), S(t) can 
be described by 

S(t) = 1 + (1/2)(1 + a) cos(2^ii)e~ 7 * 

+ (1/2)(1 - a) cos(2fi 2 0e- (n2/Ol)7t (4) 

where 7 allows for decay of the signal . The local max- 
imum at t = 2.4 fis on Fig. 0b is the 2tt : tt point at which 
ion 1 has undergone a 27r-pulse while ion 2 has undergone 
a 7r-pulse resulting in | ||)|0) — > | |T)|0)- Driving a ir : ir 
pulse on the co-propagating carrier transforms | 4t)|0) to 
I Tl)|0) and I li)\0) to | TT)|0), completing generation of 
all four internal basis states of Fig. |l|b. 

Now consider the levels coupled by the first rsb shown 
in Fig. gb. If we start in the state |^(0)) = | jT)|0) and 
drive on the (stretch mode) rsb for time t, 
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where G = (0|_ + Vl\_) 1 / 2 with Cl^ from Eq. | with 

n = 0. The phase factor <j> = Sk ■ (xi — x 2 ) depends on 
the spatial separation of the ions and arises because each 
ion sees a different phase in the x travelling-wave Raman 
field. The ion- ion spacing varies by SI « 100 nm over 
the range of Uq cited above (<^ = for Uq = 16.3 V and 
<j) = tt for Uq = 12.6 V, with d(f>/dUo in good agreement 
with theory). For Gt — tt, the final state is ijje(^>) from 
Eq. |[ Note that f2j = (y/2 + 1)0,2 would generate the 
Bell states (but we would not have access to the initial 
state I |t); since are fixed throughout an experiment). 

We now describe our two-ion state-detection proce- 
dure. We first prepare a two- ion basis state \kl), ap- 
ply the detection beam D2 for a time w 500/zs and 
record the number of photons m detected in time Td- We 
repeat this sequence for TV « 10 4 trials and build a his- 
togram of the photons collected (Fig. ||). To determine 
the populations of an unknown state, we fit its histogram 
to a weighted sum of the four reference histograms with 
a simple linear least-squares procedure. 

We observe that the | ||) count distribution (Fig. 
^a) is not a single peak at m = 0, corresponding to 
the expected zero scattered photons. Counts at m = 1 
and in = 2 are due to a background of 200-400 pho- 
tons per second. The counts in bins m > 2 (which 
account for ~ 10% of the area) are due to a depump- 
ing process in which D2 off-resonantly drives an ion out 
of I I) ultimately trapping it in the cycling transition. 
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We approximately double the depumping time by ap- 
plying two additional Raman "shelving" pulses (| |) — > 
2 S 1/2 \F = 2,m F = 0) -^ 2 S 1/2 \F = l,m F = -1); | |) 
unaffected) after every state preparation. Nevertheless, 
this results in an average difference of only 10-15 detected 
photons between an initial | |) and | f) state, as shown 
in Fig. H [p3| . The distributions associated with | |f) 
,| 11) and | ||) are non-Poissonian due to detection laser 
intensity and frequency fluctuations, the depumping de- 
scribed above and | |) — > | f) transitions from imperfect 
polarization of D2. 

One may ask: what is our overall two-ion state- 
detection efficiency on a per experiment basis? To ad- 
dress this issue, we distinguish three cases: 1) | ff)j 2) 
I Tl) or I IT); 3) | ||). Now define case 1 to be true 
when to < 3, case 2 when 3 < m < 17, and case 3 when 
to > 17. This gives an optimal 80% probability that the 
inferred case (1, 2, or 3) from a measured to in a single 
experiment is the actual case. 

We have generated states described by density opera- 
tors p^ in which the populations (diagonals of p^) are 
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measured to be Pif 

0.05. To establish coherence, consider 
singlet state ipg which has P^ 



0.4, P Ll « 0.15, P 
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TT ~ 
first the Bell 

= 1/2. Since 



tpg has total spin J — 0, any J-preserving transfor- 
mation, such as an an equal rotation on both spins, 
must leave this state unchanged, whereas such a rota- 
tion on a mixed state with populations P|f = Pfj. = 1/2 
and no coherences will evolve quite differently. We ro- 
tate both spins trough an angle 9 by driving on the 
co-propagating carrier for a time t such that 9 = Cl c t. 
Fig. ^a shows the time evolution of an experimental 
state which approximates the singlet Bell state. Con- 
trast this with the approximate "triplet" state shown 
in Fig. ^3. More quantitatively, the data show that 
p ± is decomposed as — C\ip B )(ip B \ + (1 — C)p m in 
which p rn has no coherences which contribute to the mea- 
sured signal (off diagonal elements connecting | ||) with 
| |T) and | TT) with | ||)^, and C = 0.6 is the con- 



trast of the curves in Fig. m This leads to a fidelity of 



(Vl|p ± |^) = (nT + ^U + C)/2«0.7. 

The non-unit fidelity of our states arises from several 
technical factors. The first is Raman laser intensity noise 
which gives rise to a noise-modulated Rabi-frcquency. 
The second is a second-order (in -q) effect on f2j due to 
the motional state of the COM mode 0J, which is not in 
the ground state at all times pb]] . These effects can be 
seen in Fig. |2|b as a decay envelope on the data (mod- 
eled by 7 of Eq. 0) a nd cause a 10% infidelity in initial 
state preparation P4| . This initial imperfection in state 
preparation, the contribution of the above factors on the 
rsb pulse and a first order effect due to imperfect ground- 
state preparation of the stretch mode are responsible for 
the rest of the infidelity. 

The micro-motion-induced selection of Rabi frequen- 
cies as here demonstrated is sufficient to implement uni- 
versal quantum logic with individual addressing la]. To 



isolate ion 1, we arrange the trap strength and static 
electric field so that ion 1 is at the rf null position 
(Oi = tt c Jo(0) — f2 c ) and ion 2 is at a position such 

that f^2 = Q c Jo(\8k\£,2) — 0. To isolate ion 2, we drive 
on the first motional sideband of the rf-micro-motion by 
adding S1t/27t = ±238 MHz to the difference frequency 
of the Raman beams resulting in fix = ^c^i(O) = 
and f^2 = Q c Ji(\8k\£,2) — 0.519n c . This provides a 
means of individual addressing, with which the Cirac- 
Zoller scheme can be implemented for two ions. 

In conclusion, we have taken a first step in the quan- 
tum preparation and manipulation of entangled states 
of multiple trapped ions — a step which is crucial for 
quantum computations with trapped ions. We have en- 
gineered entangled states deterministically, that is, there 
is no inherent probabilistic nature to our quantum entan- 
gling source. We have developed a two-ion state-sensitive 
detection technique which allows us to measure the diag- 
onal elements of the density matrix p^ of our states, and 
have performed transformations which directly measure 
the relevant off-diagonal coherences of p^ . 
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FIG. 2. (a) Normalized i-carrier Rabi frequencies f2i/Sl c of 
each of two ions as a function of center-of-mass displacement 
from the rf-null position d. The solid curves are Eq. ^ where 
the distance between the maxima of the two curves sets the 
scale of the ordinate, based on the known ion-ion spacing 
of I » 2.2/im at u x /2tt = 8.8 MHz. (b) Example of Rabi 
oscillations with f2i = 2D.2 . A fit to Eq. ^ determines that 
D. 1 /2ty = 2Q 2 /2ty « 225 kHz, 7/2tt « 6 kHz and a « -0.05. 
The arrow in (a) indicates the conditions of (b). 




FIG. 3. Photon-number distributions for the four basis 
qubit states. Plotted in each graph is the probability of oc- 
currence P(m) of m photons detected in 500 jj,s vs. m, taken 
over ~ 10 4 trials. Note the different scales for each graph. 
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FIG. 1. (a) Relevant 9 Be + energy levels. All optical 
transitions are near A = 313 nm, A/2ir = 40 GHz and 
LOo/2n = 1.25 GHz. Rl-3: Raman beams. Dl-3: Doppler 
cooling, optical pumping and detection beams, (b) The inter- 
nal basis qubit states of two spins shown with the vibrational 
levels connected on the red motional sideband. The labeled 
atomic states are as in (a); n is the motional-state quantum 
number (note that the motional mode frequency (jj str <Si uio). 
£li± axe the Rabi frequencies connecting the states indicated. 
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FIG. 4. Probabilities Pjf + -Pu an d -Pu + At as a function 
of time i driving on the co-propagating carrier, starting from 

(a) the "singlet" V> e (0) and (b) the "triplet" ipe.{it) entangled 
states. The equivalent rotation angle is 2Q. c t (fi c /27r ~ 200 
kHz for these data). The solid and dashed lines in (a) and 

(b) are sinusoidal fits to the data, from which the contrast is 
extracted. 
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